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Abstract. Resonance states of a two-electron quantum dot are studied using a 
variational expansion with both real basis-set functions and complex scaling methods. 
We present numerical evidence about the critical behavior of the density of states in 
the region where there are resonances. The critical behavior is signaled by a strong 
dependence of some features of the density of states with the basis-set size used to 
calculate it. The resonance energy and lifetime are obtained using the scaling properties 
of the density of states 
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1. Introduction 

Resonance states are slowly decaying scattering states characterized by a large but 
finite lifetime [H [2] . There is a host of established methods that allow the calculation 
of the energy and lifetime associated to the resonance p], EJ [3], |4j [5j El [7] • In many 
cases of interest where complex scaling (analytic dilatation) techniques can be applied, 
resonances energies and lifetimes show up as real and imaginary part, respectively, 
of isolated complex eigenvalues of the rotated Hamiltonian [TJ. Though conceptually 
simple, complex scaling is difficult to implement in molecular systems, for this kind 
of systems the analytical continuation is obtained using absorbing potential methods 
[8], which can be roughly grouped in two categories: the absorbing potentials based 
on negative imaginary potentials (NIPs), and the potentials derived from the complex 
scaling based theories (see [8] and References therein). 

Despite the wealth of methods available to deal with the calculation of the resonance 
energy the subject constantly receives attention and new methods are proposed. The 
new proposals arise, mainly, because almost all the methods have some drawbacks or 
because new concepts are applied to the problem. 

Recently Pont et al [9] and Ferron et al [10] have used the Fidelity [UJ and the von 
Neumann entropy [12], respectively, to obtain the real part of the resonance energy of 
a two electron quantum dot. Both concepts, borrowed from the quantum information 
theory, provided accurate methods for the calculation of resonance-state properties. In 
particular the method proposed by Pont et al [9] have some points in common with 
the work of Kar and Ho [7J , both methods rest on the availability of a monoparametric 
family of basis sets, where each basis set corresponds to a particular value of a free 
parameter. 

Kar and Ho [7J calculate a variational approximation for the problem which provides 
an approximate spectrum. With the spectrum they construct the density of states 
(DOS) of the problem and, after this step, it is possible to find the complex eigenenergy 
of the resonance fitting the DOS. Their method is quite general but the fit that is 
necessary to obtain the resonance energy is by no means unambiguously defined. In 
contradistinction the methods proposed by Pont et al [9] do not require any fitting, but 
to produce the resonance energy for a given problem (with its corresponding parameters) 
it is necessary to find an adequate value of the basis-set free parameter by bisection. 

In this work we will present a method to find the energy of a resonance state using 
the DOS but without resort to any fitting or bisection. We will study the resonance 
state that arise in a two electron quantum dot when its ground state looses stability and 
enters into the continuum. In Section |2] we briefly present the model and the variational 
method used to find an approximate spectrum. In Section[3]the DOS is obtained and the 
energy and lifetime of the resonance state are founded using finite size scaling techniques. 
Finally, in Section [4] we discuss our results and draw some conclusions. 
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2. The model and basic results 

There are many models of quantum dots, with different symmetries and interactions 
(see [U [10] and References therein). In this work we consider a model with spherical 
symmetry, with two electrons interacting via the Coulomb repulsion. The main results 
should not be affected by the particular potential choice as it is already known that the 
near threshold behavior and other critical quantities (such as the critical exponents of 
the energy and other observables) are mostly determined by the range of the involved 
potentials [13]. Therefore to model the dot potential we use a short-range potential 
suitable to apply the complex scaling method. After this considerations we propose the 
following Hamiltonian H for the system 

H = -fvl - ^Vl 2 + V(n) + V(r 2 ) + y*— (1) 
2m 2m |r 2 — ri| 

where V(r) = — (Vo/tq) exp (— r/r ), rj the position operator of electron i = 1,2; r 
and Vo determine the range and depth of the dot potential. After re-scaling with tq, in 
atomic units, the Hamiltonian of Eq. (P) can be written as 

H = -\vl, - ~V* 2 - V e^ - V e-r> + j-A— ( 2 ) 

where A = r . 

We choose the exponential binding potential to take advantage of its analytical 
properties. In particular this potential is well behaved and the energy of the resonance 
states can be calculated using complex scaling methods. So, besides its simplicity, the 
exponential potential allows us to obtain independently the energy of the resonance 
state and a check for our results. The threshold energy, e, of Hamiltonian Eq. that 
is, the one-body ground state energy can be calculated exactly [U] and is given by 

■J2V2S (VWo) = 0, (3) 

where J v {x) is the Bessel function. 

The discrete spectrum and the resonance states of the model given by Eq. (j2J) can 
be obtained approximately using C? variational functions prj, [6]. So, if |^(1,2)) are 
the exact eigenf unctions of the Hamiltonian, we look for variational approximations 

M 

|^(1,2)) ~ 1^.(1,2)) = J2 C ^\^) » c ? ) = (c (j) )* ; J = 1,---,M.(4) 

i=l 

where the |$,) must be chosen adequately and M is the basis set size. 

The functions are chosen as eigenfunctions of the total angular momentum with 
zero eigenvalue. The radial part of the basis functions are constructed by symmetrization 
of the one-body Slater functions 



The scaling of the density of states 



4 



With these constraints the functions have only three relevant quantum numbers 
ni, n 2 and I, where ni and n 2 are related to the radial part of the function, and I to the 
angular one [9]. The basis set is characterized by a non-linear parameter a, then we use 
the notation — > \ni, n 2 ; I; a). So, we get that 

^(l,2))=^c«£ ) |n 1 ,n 2 ;/;a), (6) 
where ni > n 2 > / > 0, then the basis set size is given by 

N n\ ri2 , 

M = E EE 1 = ^(N+l)(N + 2)(N + 3), (7) 

ni=0n 2 =0 2=0 

so we refer to the basis set size using both N and M. The matrix elements of the kinetic 
energy, the Coulombic repulsion between the electrons and other mathematical details 
involving the functions \ni,n 2 ;l;a) are given in references [9], [15], [16] . 
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Figure 1. (color on-line) (a) the figure shows the behavior of the variational 
eigenvalues Ej(A) (black lines) for N = 14 and non-linear parameter a = 2. The red 
dashed line corresponds to the threshold energy e ~ —1.091. Note that the avoided 
crossings between the variational eigenvalues are fairly visible, (b) The figure shows 
the same variational eigenvalues that (a) (black lines) and the energy calculated using 
the complex scaling method (green line) for a parameter 9 = n/10. 

Resonance states have isolated complex eigenvalues, E res = E r — iT/2, T > 0, 
whose eigenfunctions are not square-integrable. These states are considered as quasi- 
bound states of energy E r and inverse life time T. For the Hamiltonian Eq. (J2J), the 
resonance energies belong to the interval (e,0) [2]. 

The resonance states can be analyzed using the spectrum obtained with a basis 
of C 2 functions (see [10] and References therein). The levels above the threshold have 
several avoided crossings that "surround" the real part of the energy of the resonance 
state. The presence of a resonance can be made evident looking at the eigenvalues 
obtained numerically. Figure [1] shows a typical spectrum obtained from the variational 
method. This figure shows the behavior of the variational eigenvalues as functions 
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of the parameter A, the results were obtained using N = 14 and a = 2.0. The value of a 
was chosen in order to obtain the best approximation for the energy of the ground state 
in the region of A where it exists. The figure shows clearly that for A < Xth ~ 1.54 there 
is only one bound state. Above the threshold there is not a clear cut criteria to choose 
the value of the non-linear parameter, and the variational approximation provides a 
finite number of solutions with energy below zero. However, it is possible to calculate 
E r (\) calculating for many different values of the variational parameter (see Kar 
and Ho 0). 

Figure |2] shows the numerical results for the second eigenvalue E^\\)^ for different 
values of the variational parameter a. As can be seen from Figure [2} there are three 
regions, in each one of them the curve for a given value of a is basically a straight line 
and the slope is different in each region. A feature that appears rather clearly is that, 
for fixed A , the density of levels per unit energy is not uniform, despite that the curves 
E^ l \\) are drawn for forty equally spaced ccj's between a = 2.0 and a = 6.0. This 
fact has been observed previously pT] and the DOS can be written in terms of two 
contributions, a localized one and an extended one. The localized DOS is attributed to 
the presence of the resonance state, conversely the extended DOS is attributed to the 
continuum of states between (e, 0). 




Figure 2. (color on-line) The second variational state energy vs A, for different values 
of the variational parameter a. From bottom to top a increases its value from a = 2 
(dashed blue line) to a = 6 (dashed orange line). The real part of the resonance 
eigenvalue obtained using complex scaling ( 9 = w/10) is also shown (green line). 
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3. The density of states 

The localized DOS p(E) can be expressed as [T7j 



P(E) 



dE(a) 



-i 



da 

Since we are dealing with a variational approximation, we calculate 



a i+ i - cti-x 



(8) 



(9) 



In complex scaling methods the Hamiltonian is dilated by a complex factor i] = ae~ l9 . 
As was pointed out long time ago by Moiseyev and coworkers, the parameter a is 
equivalent to our parameter a [18]. Besides, the DOS attains its maximum at optimal 
values of a and E r that could be obtained with a self-adjoint Hamiltonian without 
using complex scaling methods [19]. So, locating the position of the resonance using 
the maximum of the DOS is equivalent to the stabilization criterium used in complex 
dilation methods that requires the approximate fulfillment of the complex virial theorem 

Figure [3] shows the typical behavior of pj(E) = p(Ej at \X)) for several eigenenergies, 
with A = 2.25. The real and imaginary parts of the resonance's energy, E r (X) and 
Ei(X) = — T/2 respectively, can be obtained from p{E), see for example [7] and references 
therein. 

The values of E r (X) and T(A) are obtained performing a nonlinear fitting of p{E), 
with a Lorentzian function, 

p(E) = Po + 7 [(^-^ /2 +(r/2)r (10) 

One of the drawbacks of this method results evident: for each A there are several 
Pj{E) (in fact one for each variational level, see Figure [3]), and since each Pj{E) provides 
a value for E 3 r (X) and T-^A) one has to choose which one is the best. Kar and Ho [7] solve 
this problem fitting all the Pj{E) and keeping as the best values for E r (X) and T(A) the 
fitting parameters with the smaller \ 2 value. At least for their data the best fitting (the 
smaller \ 2 ) usually corresponds to the larger n. This fact has a clear interpretation, 
if the numerical method approximates E r (X) with E^\X) a large n means that the 
numerical method is able to provide a large number of approximate levels, and so the 
continuum of states between (e, 0) is "better" approximated. 

In most numerical studies it is assumed that using the largest possible basis set size 
leads to the best results. Indeed, most usually this assumption is correct but taking in 
consideration only the results obtained for the largest basis set size does not provide 
any insight about the accuracy of the result. 

The Finite Size Scaling method (FSS) provides a systematic way to analyze the 
data generated using different basis set sizes and has been extensively used in atomic 
physics [21], statistical mechanics [22], etc. The FSS method requires a scaling function 
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Figure 3. (color on-line) The DOS p(E) for A = 2.25 and basis set size N = 14. The 
results were obtained using Eq. (|8]) and correspond, from top to bottom, to the second 
(black line), third (dashed red line), fourth (green line) and fifth (dashed blue line) 
levels. 

that reduces the data corresponding to different basis set sizes to a single "universal" 
set. Moreover, as we will show, the FSS allows to obtain the real and imaginary parts 
of the resonance energy using any Pj(E). 

We chose p%{E) to apply the FSS method, but using any other pj should lead to the 
same conclusions. As we have said previously the functional form of the scaling function 
is unique, independently of the actual value of A. Figure H](a) shows p${E) for different 
values of A and several basis set sizes. For a fixed value of N, say N = 16 the Figure 
11(a) shows seven peaks, each one corresponding to a given value of A. The values of 
A were chosen to sample different regimes of the resonance, from small values -sharp 
resonances- to larger ones -broad resonances. The broadening of the peaks corresponds 
roughly to the broadening of the resonance or, equivalently, to larger values of T. 

Figure H](b) shows the scaled data of panel (a), for each A the curves corresponding 
to different values of N collapse into a single one which is independent of N. 

The form of the scaling function can be inferred from the behavior of the maximum 
value of the function pz{E) for different basis set sizes. Figure [5] shows the behavior 
of ln(p 3max ) vs ln(iV), for N large enough the maximum values behave as a power law. 
With this numerical evidence we propose 

p{E) — ► N-^p(E) (11) 

as the scaling function for the DOS. As can be seen in Figure H] (b), the collapse of 
the data is excellent. The exponent /3 depends on A and its value must be found 
for each case. So, for a given A the real part of the complex resonance energy 
corresponds to the value where the DOS attains its maximum, E r (X) = E Pmax , where 





Figure 4. (color on-line) (a) The DOS pz{E) for (from left to right) A = 
2.25, 2.50, 2.75, 3.00, 3.25 and 3.50 and different basis set sizes. The results were 
obtained using Eq. (J8j) - (b) Data collapse obtained using the scaling function 



Pmax = max£ p(E) = p(E Pmax ). At least from our data it is clear that the localization 
of the maximum of the DOS is a very slowly function of the basis set size, see Figure 
(b). 
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Figure 5. (color on-line) (a) log — log plot of the maximum of the density of states 
against the basis-set size N. (b) E r vs. 1/N. For both panels the data were obtained 
with A = 2.25. 



On the following we will argue about the form of the scaling function necessary to 
obtain the imaginary part of the resonance energy, T/2. As we have said previously, the 
real part of the resonance energy is located between (e, 0), and as log as T — > when 
A — > A+, the scaling function proposed for the imaginary part of the resonance energy 
should be zero at A c or, equivalently, E r (X c ) = e, and Ei(\ c ) = — r(A c )/2 = 0. Besides, 
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in some examples the life time of a resonance is proportional to the DOS, r oc -, so we 
tested our numerical results using 



E i l N \\) = -K (EW{\)-E C )W 



AN) ■ 

ymax 



(12) 



where we have made explicit the superscript N to emphasize the dependency with the 
basis set size, kappa is a positive constant, and 5 is a exponent to be determined. As 
a first guess we used S = 0.8, since this figure is very close to the exponents founded 
studying the scaling of the maximum of the DOS. Then, we observed that 



£f°(A") E^ a> {\" 



(CS), 



(13) 



i. e. the ratio between the imaginary parts of two resonance energies, corresponding to 
A" and A' chosen arbitrarily, calculated using the Complex Scaling (CS) method and 
Eq. (Tl2l) is fairly the same. Regrettably this finding does not provide the value of 
the constant k. Just to compare with the result of Complex Scaling methods we 



calculated n imposing that E\ N \\) = E^ S, {X) for one particular A. Figure |6] shows 
the values obtained for e\ obtained using the recipe described above and Eq. ( TT2l) . 
the agreement between values corresponding to different basis set sizes is excellent, and 
the same can be said with respect to the values obtained from Complex Scaling, see 
Figure [6] (a). Despite that the value of k was obtained using a particular value of N and 
A our numerical findings strongly support the functional form proposed in Eq. ffl2l . in 
Figure |6] (b) we present a comparison between the E^ N \\) y s and the Complex Scaling 
results, the agreement is better than 1% for all the values of A and N analyzed. 
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Figure 6. (color on-line) (a) Resonance pole string curve obtained from complex 
scaling and basis-set size scaling. A mean value of S = 0.8 was chosen for the scaling 
exponent, (b) The energy discrepancies between the two methods. The relative error 
is less than 1% for the calculated values. 
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4. Summary and conclusions 

So far we have presented numerical evidence of the critical behavior of the density of 
states near the peak associated to the presence of a resonance state. The critical behavior 
should be understood in the sense of FSS method, i.e there is a strong dependency of 
the density of states with the basis set size used to obtain it. The critical behavior 
also manifest itself through the scaling properties of some quantities, in particular the 
maximum value attained by the density of states for fixed A and N. 

Despite that we have been able to propose a scaling function that provides the 
imaginary part of the resonance energy, the proposal is made only in phenomenological 
terms and, more importantly, relays in other methods to be quantitative. Our efforts 
will be directed to put the scaling law Eq. (fl2|) over a more formal basis, which hopefully 
will allow us to obtain the constants k and 5 using only scaling arguments. 

The critical behavior of the density of states studied in this work, together with 
the behavior observed in the fidelity [9], points that the analogy between the "line" of 
resonance states and a "critical line" (in the sense of critical phenomena) could be more 
exploited and it is worth of further exploration. 
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